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ABSTRACT 

Although the occurrence of steady supercritical disk accretion onto a black hole has been speculated 
about since the 1970s, it has not been accurately verified so far. For the first time, we previously 
demonstrated it through two-dimensional, long-term radiation-hydrodynamic simulations. To clarify 
why this accretion is possible, we quantitatively investigate the dynamics of a simulated supercritical 
accretion flow with a mass accretion rate of ^ 10^Le/c^ (with Le and c being, respectively, the 
Eddington luminosity and the speed of light). We confirm two important mechanisms underlying 
supercritical disk accretion flow, as previously claimed, one of which is the radiation anisotropy arising 
from the anisotropic density distribution of very optically thick material. We qualitatively show that 
despite a very large radiation energy density, £'o^lO^-^E/47rr^c (with r being the distance from the 

black hole), the radiative flux Fq ^ cEq/t could be small due to a large optical depth, typically 
T 10'^, in the disk. Another mechanism is photon trapping, quantified by vEq, where v is the flow 
velocity. With a large and Eq, this term significantly reduces the radiative flux and even makes 
it negative (inward) at r < TOrg, where is the Schwarzschild radius. Due to the combination of 
these effects, the radiative force in the direction along the disk plane is largely attenuated so that the 
gravitational force barely exceeds the sum of the radiative force and the centrifugal force. As a result, 
matter can slowly fall onto the central black hole mainly along the disk plane with velocity much less 
than the free-fall velocity, even though the disk luminosity exceeds the Eddington luminosity. Along 
the disk rotation axis, in contrast, the strong radiative force drives strong gas outflows. 
Subject headings: accretion: accretion disks — black hole physics — hydrodynamics — radiative 
transfer 



1. INTRODUCTION 

Recently, very bright objects which may be under- 
going supercritical (or super-Eddington) accretion flows 
have successively been found. Good e xamples are ul- 
traluminous X-ray s o urces fULXs: Watarai et all lMll: 
lEbisawa et all 120031 : lOkajima et al ll2006( l. These are 
pointlike off-center X-ray sources whose X-ray luminos- 
ity significant ly exceeds the E ddington luminosity of a 
neutron star ( Fabbianol Il989f ) . Because of substantial 
variations, it is reasonable to assume that the ULXs 
are single compact obi ects powered by accretion flows 
(|Makishima et al1l2000[ ). If so, there are two possibili- 
ties to account for large luminosities exceeding the Ed- 
dington luminosity for a mass of 100 Mq: sub-critical 
accretion onto an intermediate-mass black hole (IMBH) 
and supercritical accretion onto a stellar-mass black hole. 
We support the latter possibility, since through the fit- 
ting to several XMM-Newton EPIC data of ULXs, which 
have been claimed as good IMBH candidates, we have 
found evidence of supercritical flows (Vierdavant i et all 
I2OO6I ). Another interesting group is narrow- line Seyfert 
1 galaxies fsee iBolleil |2004 for a review). Because of 
their relatively small black hole masses, they have in 
general large Eddington ratios (L/L-e with L and Le 
being the luminosity and the Eddington luminosity, re- 
spectivel y) , and some of them seem to fall in the s lim-disk 
regmies (Mineshige et al.ll2000l : lKawaguchill200a see also 



iWang et al.l[l999l ) 

Despite growing evidence indicating the existence of 
supercritical accretion flows in the universe, theoreti- 
cal understanding is far from being complete. It is well 
known that any spherically accreting object, irrespective 
of the nature of the central source, cannot emit above 
the Eddington luminosity, since otherwise significant ra- 
diative force will prevent accretion of the gas. If we ex- 
amine detailed radiation-matter interactions in the in- 
terior, however, we notice that the situation is not so 
simple. Actually, radiation produced at the very cen- 
ter cannot immediately reach the surface (i.e., photo- 
sphere), since photons generated at the center should 
suffer numerous scatterings with accreting material and 
thus take a long time to reach the surface. If the matter 
continuously falls, and if the mass accretion timescale is 
shorter than the mean travel time for photons to reach 
the surface (the diffusion timescale), photons at the core 
may not be able to go out. This is the so-called photon- 



trapp ing effect (e.g.. lBegelman|[T978tlHouck &: Chevalieil 
Il991f) . Here we define the trapping radius, inside which 
radiation-matter interaction is so frequent that photons 
are trapped within the accretion flow. Inside this trap- 
ping radius, therefore, radiative flux can be negative (i.e., 
inward). Thus, the apparent luminosity is reduced as 
compared with the case without photon trapping. Nev- 
ertheless, the concept of the Eddington luminosity is still 
valid, since far outside the trapping radius radiative flux 
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should be outward and its absolute value should be less 
than L^/A-KT^, where r is the distance from the cen- 
tral black hole, in the quasi-steady state. Radiation- 
hydrodynamic (RHD) simulations of spherically symmet- 
ric supercritical accre tion flows have been performed by 
iBurger fc KaX^ ([l983h . 

The situation may differ in the case of disk accretion, 
since the radiation field is not isotropic due to inhomo- 
geneous matter distribution. That is, matter can fall 
predominantly along the disk plane, whereas radiation 
can go out along its rotating axis, where matter den- 
sity is low. In other words, the main directions of the 
inward matter flow and the outward radiative flux are 
not parallel to each other in disk accretion, leading to 
a situation in which the radiative force does not com- 
pletely counteract the gravitational force. There is room 
for the possibility of a supercritical accretion flow with 
super-Eddington luminosity. 

Based on such an argument, many researchers have 
speculated about the occurren ce of supercritica l flows 
in disk accretion systems. Shakura fc SunvaevI ()1973f ) 
discussed the possibility of supercritical disk accretion 
based on a one-dim ensional steady model (see also 
iMaraschi et al.l I1976D . They mentioned that the mass 
accretion rate would not be steady but oscillate if the 
mass accretion rate exceeded the critical rate; other- 
wise, a part of the accreting matter might be ejected 
from the disk as a disk wind. Radiatively drive n out- 
flows fr om supercritica l disks were inves t igated bv lMeiej 
([i97i),[j ones fc RBind (p79h . andlMi (fl98l . 

Despite a long history in the study of supercritical disk 
accretion flows, the occurrence of steady supercritical 
disk accretion has not yet been accurately verifled. Simi- 
lar simulations have been performed since the 1980s, but 
all of them calculated only the initial transient phase. 
Their conclusions then were not general, since the back- 
reactions (i.e., enhanced radiation pressure), which may 
inhibit steady flow, were not accurately evaluated. Al- 
though the radiative force predominantly has an effect in 
the vertical direction, it should also have an effect on the 
material within the disk. Hence, in the direction paral- 
lel to the disk plane, the situation may be similar to or 
more severe than the case of spherical accretion. This is 
because the radiative force, together with the centrifu- 
gal force, may possibly overcome the gravitational force. 
To study the possibility of supercritical disk accretion 
flows, precise and quantitative research treating both su- 
percritical disks and outflows, and which also takes into 
c onsideration niulti-di mensional effects, is needed. 

lOhsuga et al.l ()2005L hereafter Paper I) have confirmed 
the occurrence of quasi-steady supercritical disk accre- 
tion onto a black hole by two-dimensional RHD simula- 
tions. The motivation of the present study is to investi- 
gate the physical mechanisms which make supercritical 
disk accretion possible. For this purpose we examine 
quantitatively the flow motion and force fields via the 
radiative flux, rotation, and gravity of a supercritical 
disk accretion flow onto a central black hole, based on 
the two-dimensional RHD simulation data from Paper 
I. Through detailed inspection of the results, it will be 
possible to clarify the physics behind supercritical disk 
accretion flows. In §2 we plot the spatial distributions 
of several key quantities which control flow dynamics. A 
discussion is given in §3. 



2. DYNAMICS OF A SUPERCRITICAL DISK ACCRETION 

FLOW 

2.1. Basic Considerations 

Here we adopt the spherical coordinates (r, 0, (/?), with 
9 — it/ 2 corresponding to the disk plane and the ori- 
gin being set at the central black hole. The momentum 
equation for the radial component of the flow is 

"J^ — /grav + /rad + /prcs + /cent' (1) 

Here Vr is the radial component of the flow velocity, 
d/dt = d/dt + Vrd/dr is the Lagrangian derivative, 
/grav is the gravitational force by the central black hole, 
/Jad — xFq/c is the radiative force with x being ab- 
sorption and scattering coefficients, c being the speed of 
light, and Fq representing the radial component of the 
radiative flux in the comoving frame, /pros is the radial 
component of the gas pressure force, and /cent — 'v'^/i^ 
is the radial component of the centrifugal force. We 
adopt a pseudo-Newtonian potential, thus having /grav = 
—GM/ (r — rg)^ wit h M being the mass of the black hole 
and rs = 2GM/c^ (iPaczvhskv fc WiitalflMih . 

To proceed, it is important to distinguish two different 
views: the view from an observer comoving with the ac- 
creting gas (i.e., the comoving frame) and the view from 
an observer standing at infinity (the rest frame). The ith 
component of radiative fiux in the comoving frame, Fq, 
is related to that in the rest frame, F'', on the order of 
v/c, 

F' ^F^+v,Eo + v,P;\ (2) 

where Eq and Pq'' are the radiation energy density and 
the radiation stress tensor in the comoving frame, re- 
spectively. Here the final term in equation ([2]) is equal 
to ViEo/S in the optically thick diffusion limit. 

In §2.3 we show that the gas accretes toward the black 
hole in the disk region. In this region, even if the radia- 
tive flux in the comoving frame is positive, Fq > (i.e., 
outward flux), F^ can be negative, F'' < (inward flux), 
because Vr < 0. The term VrE^ thus represents photon 
trapping. In §2.5 we explicitly show this effect in the 
simulation data. 

In contrast, the ffow dynamics in the direction perpen- 
dicular to the disk plane is distinct. In particular, the 
centrifugal force has no effect in this direction. In addi- 
tion, the matter density is significantly smaller along the 
rotation axis. Hence, the radiative flux can be much more 
effective in the vertical direction, thus driving strong out- 
flows, as we see in §2.3. 

2.2. Overview of the Simulated Flow 

Our present analysis is based on the two-dimensional 
RHD simulation data from Paper I. These are the first 
simulations of supercritical accretion flow in the quasi- 
steady state. Although the research histor y of such sim- 
ulations stems back to the late 1980s, when lEggum et ahl 
f(1987.i ) performed numerical simulations for the first time, 
their cal culations were restricted to the first few seconds 
(see also'Klevl ll98l lOkuda et al.lfl997t iKlev fc LinlfT999l : 
lOkuda 200^ Back-reactions were not fully taken into 
account in their simulations. In our simulations pre- 
sented in Paper I, matter is added continuously from 
the outer disk boundary at r = 500 rs to the initially 
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Fig. 1. — Color contours of matter density (left) and radiation c 
these values are time-averaged over t = 20 — 50 s. 

(nearly) empty space at the rate of lO'^L^/c^. The in- 
jected mass has a small specific angular momentum so 
that it first free falls but then forms a rotating disk at 
r ~ lOOrg. Mirror symmetry is assumed with respect 
to the disk (equatorial) plane. The a-prescription of the 
viscosity is employed, in which the viscosity is set to be 
proportional to the total pressure in the optically thick 
limit and the gas pressure only in the optically thin re- 
gion. Hence, the viscosity is more effective in the disk 
region than in the outflow regions above and below the 
disk. The viscosity parameter is assumed to be a = 0.1. 
Except for the radial-azimuthal component, all compo- 
nents of the viscous stress tensor are set to be zero. 

The radiative transfer was solved un- 
der the flux-limited diff usion approximation 
(jLevermore Sz Pomranin j Il981l ). The central ob- 
ject is taken to be a non rotating stellar-mass black 
hole (M = IOMq), generating a pseudo-Newtonian 
potential. Only 10% of the inflowing material finally 
reaches the inner boundary (Srg); i.e., 90% of the mass 
input gets stuck in the dense, disk like structure around 
the equatorial plane, or transforms into the known 
collimated high-velocity outflows perpendicular to the 
equatorial plane or into low-velocity outflows with wider 
opening angles. The resulting luminosity is about 3 
times larger than the Eddington luminosity. 

Figure [T] indicates the time-averaged contours of the 
matter density {left panel) and the radiation energy den- 
sity {right panel) on the R-z plane {R = rsinO and 
z = rcos9), where both axes are normalized by the 
Schwarzschild radius {rs). Note that these panels are 
similar to Figures 3 and 5 of Paper I, except that Fig- 
ure [1] of the present paper has contours which are time- 
averaged over t = 20 — 50 s, whereas those in Paper I 
are snapshots at t — 10 s (i.e., in a quasi-steady phase). 
The high concentration of gas and radiation in the region 
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gy density (right) distributions as a function of R and z. Note that 

near the black hole on the equatorial plane is clear. 

2.3. Dynamics of the Simulated Flow 

First we plot the time-averaged radiation energy den- 
sity normalized by L-e/ Anr'^c in Figure [2l We find that 
this value is larger than unity in most regions. In partic- 
ular, we have £^o^l0^iE/47rr^c in the disk region. This 
means that there exists so large a number of photons that 
if we simply assume Fq ~ cEq as in the optically thin ve- 
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Fig. 2. — Same as Fig. 1, but with the radiation energy density 
normalized by LE/47rr^c. 
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Fig. 3. — Same as Fig. 1, but with the color contours of the ratio of the radial component of the radiative force to the gravitational force 
(left) and that of the sum of the radiative force and the centrifugal force to the gravitational force (right). 

gion, the radiative force should highly exceed the gravity, 
which inevitably prevents the inflow motion. The simu- 
lation data, however, do show inflow motion (discussed 
below) . 

Figure [3] illustrates the time-averaged force balance 
in the radial (r) directions with various ^-values from 
the origin. (It should be noted that since force is a 
vector, we need to specify the direction to draw such 
contours. Here we only plot the radial component of 
the forces.) Figure [3] {left) represents the ratio of the 
radial component of the radiative force to the gravita- 
tional force. Especially surprising is that most regions 
are turquoise, indicating that the radiative force is out- 
ward and less than the gravitational force. Why is the 
radiative force so weak, even though the radiation energy 
density is very large, Eq>1Q'^ L^/Airr'^c (see Fig. [2|), and 
even though the total luminosity exceeds the Eddington 
luminosity, L > Le? This is because the disk is very 
optically thick. Although the absolute value of the ra- 
diative flux is equal to cE^ in the optically thin limit, 
it is roughly given by cEq/t, and hence, Fq can be very 
small, Fo <C cEq, in the optically thick region (r ;» 1), 
where r is the optical thickness. Approximately, we find 
T = pxr ~ 10^(yo/10~''-^gcm~'^)(r/30r5), where p is the 
gas density. Thus, flux is attenuated a great deal. We 
thus have a relatively small radiative flux within the disk 
in spite of a large Eq (see §2.4 for more detail). There 
are two exceptions: the color of the region in the vicinity 
of the black hole is blue, which means that the radiative 
force is inward. In this region, a large number of pho- 
tons are swallowed by the black hole, Fq < 0, so that the 
gas is accelerated inward. The region near the rotation 
axis, in contrast, has a white color, which means that the 
radiative force exceeds the gravitational force. 

To discuss how matter flows, we also need to consider 
the centrifugal force. (Note that the gas pressure is neg- 
ligible in the present case.) Figure^ (right) fllustrates 




Fig. 4. — Same as Fig. 1, but with the radial velocity normalized 
by the free-fall velocity. 

the ratio of the sum of the radial components of the ra- 
diation and centrifugal forces to the gravitational force. 
Certainly, the repulsing force is enhanced along the disk 
plane. However, this force ratio is still around unity in 
most regions except in the very vicinity of the black hole. 
That is, the gas is nearly in balance of forces in wide 
regions. The regions near the rotation axis again have 
white colors, which means that the radiative force dom- 
inates, as in Figure [3] (lejf) (Note that the centrifugal 
force is not effective near the rotation axis.) 
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Since there are regions with a negative effective force 
in the radial direction, we expect regions with negative 
velocity (i.e., inflow). Figure [J] shows the spatial distri- 
bution of the radial component of the matter velocity in 
units of the free-fall velocity. Interestingly, the flow is 
inward in the disk region, but the infall velocity is much 
smaller than the free-fall velocity. Along the rotation 
axis, conversely, the flow is outward, and its magnitude 
is larger than the free-fall velocity. 

Let us have a more quantitative discussion based on the 
one-dimensional distributions of these quantities, since it 
is difficult to obtain their precise values from the color 
contours. We calculate several quantities along two ra- 
dial lines: one at 0.457r (close to the equatorial plane) 
and the other at 9 = 2.3 x 10~^7r (nearly along the ro- 
tation axis). They are illustrated in the top and bottom 
panels of Figure O respectively. 

We see in the top panel that the radiative force is 
not negligible, although it is smaller than the centrifugal 
force. The sum of the radiative and centrifugal forces is 
nearly balanced with the gravitational force. Radiation 
does work as a "radiation cushion," which decelerates 
the accretion of the gas in cooperation with the rotation 
(centrifugal forc e) . Deceleration via radiativ e flux was 
also reported bv lBurger fc Katd (| 19801 . 119831 ). Figure [5] 
(top) also shows that the radial velocity is much smaller 
than the free-fall velocity and only barely below zero (in- 
flow), although there also exist regions with Vr ^ (no 
inflow) . 

In contrast. Figure [5] {bottom) shows that the ratio of 
the radiative force to the gravitational force grows as 
the distance increases. The outward velocity thus grows 
rapidly upward. We confirm again that the strong out- 
flow around the rotation axis is produced by the strong 
radiative force. However, the velocity is negative in the 
very vicinity of the hole. This means that the radiation 
which produces the huge radiative force in the vertical 
direction does not originate from the region just outside 
the event horizon but from the region with Eq at maxi- 
mum [at (5 — 10)r5]. 

Here we note that the angular momentum is very small 
around the rotation axis, leading to a negligible cen- 
trifugal force (dotted line). This is because a significant 
amount of material exists even inside the radius of the 
marginally stable orbit, R < 3rs, in the supercritical 
flow, so the matter can lose angular momen tum there 
(|Kato et al.lll998HWatarai fc Miiieshigell2003f ). This sit- 
uation contrasts with that of the subcritical flow, in 
which there is a sharp density drop at i? = 3rs- 

Our computational domain covers exactly down to the 
rotation axis. Numerical simulations show that most of 
the outflow matter is blown away into an oblique direc- 
tion by the radiative force in cooperation with the cen- 
trifugal force. 

2.4. Geometric Effect 

Here we consider the reason why the radiative force 
in the radial direction does not dominate over the grav- 
itational force inside the disk, even though Eq is larger 
than lO^L^/ATTr^c as shown in Figure [51 The radial 
component of the radiative flux is evaluated as Fq ~ 
(c/3px) dEo/dr in the disk region. It is attenuated since 
the disk is dense and very optically thick (t :s> 1). That 
is, even though the radiation energy density itself is large, 




Fig. 5. — time-averaged one-dimensional profiles of the radiative 
force {solid lines), the centrifugal force {dotted lines), their sum 
{dashed lines), and the radial velocity {long- dashed lines) for the 
case with a large polar angle, i = 0.457r {top), and with a small 
angle, i = 2.3 X 10~^7r {bottom). Note that all the forces are 
normalized by the gravitational force, while the radial velocity is 
normalized by the free-fall velocity. In both panels the two thin 
horizontal lines indicate the values of and 1. 

the high matter density significantly suppresses the ra- 
diative flux inside the disk. However, the high density 
(large optical thickness) is not the only condition for the 
occurrence of supercritical accretion. In the case of a 
spherical system, the radiative flux should be L/Airr^ 
even if the matter is very optically thick. The radia- 
tive force dominates over the gravitational force when 
L > Le, preventing the accretion of matter. Hence, 
there should be another reason that supercritical disk 
accretion is realized in spite of L > Le. 

Whereas the radiation field is isotropic in the spheri- 
cal case by definition, it can be anisotropic in the disk 
case. Indeed, the simulated profile of the radiation en- 
ergy density is highly anisotropic (Fig. [l] right). Unlike 
in spherical geometry, photons can escape from the less 
dense region around the rotation axis without thrusting 
through the dense disk region. Thus, the effective radia- 
tive force is attenuated in the disk region. Here we note 
that, even in the presence of an anisotropic matter dis- 
tribution, L cannot exceed Le if the medium is optically 
thin. To sum up, supercritical disk accretion with super- 
Eddington luminosity is realized due to a large optical 
thickness and anisotropy of the radiation field. 

2.5. Photon-trapping Effects 

Photon-trapping effects assist the occurrence of super- 
critical accretion. Figure [5] (top) indicates the radiative 
fluxes in the comoving and rest frames inside the disk. As 
shown in this figure, the radial component of the radia- 
tive flux in the rest frame [dotted line) is negative in the 
region of r < TOrg, whereas that in the comoving frame 
[solid line) is positive except in the vicinity of the black 
hole. This means that the radiation energy is transported 
inward via photon trapping, VrE^ < Q (see equation [5]). 
Some part of this radiation energy is swallowed by the 
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black hole, whereby the luminosity is reduced. This effect 
attenuates the radiative force, supporting the accretion. 

The anisotropy of the radiation field is enhanced by 
photon trapping. This is another important role of the 
photon trapping. The radiation energy, which is ad- 
vected inward with the matter by the photon trapping 
but not swallowed by the black hole, is transported to 
the vertical direction and finally released from the disk 
surface. It is shown in Figure 5 (bottom), which plots 
the ^-components of the radiative fluxes in the comov- 
ing {solid line) and the rest {dotted line) frames along 
a radial line at 9 = O.Stt, that the entire region is be- 
low the photosphere. In this figure we can see that the 
radiation energy is transported toward the disk surface, 
< 0, in most of the region. We also find < Fq in 
this figure. This implies that the gas motion contributes 
to the transport of radiation energy toward the disk sur- 
face, vgEo < 0. To conclude, the photon-trapping effect 
enhances the anisotropy of the radiation field, assisting 
the occurrence of supercritical accretion. 

Here we note that the gas motion in the vertical di- 
rection seems to be driven by convection in cooperation 
with the radiative force. In fact, we find that the adia- 
batic conditions for radiation-pressure-dominated matter 
are roughly satisfied along the vertical line in the disk re- 
gion, dlnP/dlnp ~ 4/3 and dlnT/dlnP - 1/4, where 
P is the total pressure and T is the temperature. 

3. DISCUSSION 

3.1. Structure of the Supercritical Disk Accretion Flow 

We summarize our simulation results in a schematic 
picture (see Fig. [71). In this figure the radiative fiux 
in the comoving frame, Fq, is positive (outward) except 
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Fig. 6. — Top: Radial components of the time-averaged radiative 
fluxes in the comoving frame (F^; solid line) and the rest frame 
{F'^; dotted line) near the equatorial plane as a function of the 
radius. Bottom: Same as the top panel, but with the polar compo- 
nents of the radiative fluxes in the comoving frame, Fq (solid line), 
and in the rest frame, F^ (dotted line) along a line with 9 = O.Stt. 



in the vicinity of the black hole. However, the radia- 
tive flux in the rest frame, F^{^ Fq + ivrEo/S), is neg- 
ative (inward) via photon trapping, VrEo < 0, in the 
trapping region. In the vicinity of the black hole, both 
Fq and F^ are negative. Matter slowly accretes inside 
the disk, since the sum of the radiative and centrifu- 
gal forces is nearly balanced with the gravitational force. 
Here we stress again that the radiative force counteracts 
the gravitational force in spite of the trapping region be- 
cause Fq > 0. Radiatively driven high- velocity outflows 
appear above and below the disk. In the very vicinity 
of the black hole, the gas is accelerated inward by the 
radiative force and the gravitational force. 

As far as the physical quantities around the equato- 
rial plane are concerned, the simulated profiles of the 
density, temperature, and radial as well as rotational ve- 
locities rougM;5;_agree_wlth of the slim-disk 
model (|Abramowicz et aLlll988D . Such features have al- 
ready been shown in Figure 11 in Paper I. However, only 
about 10% of the injected mass can accrete onto the black 
hole, and an almost equal amount of matter is ejected as 
high- velocity outfiows. The mass accretion rate is not 
constant in the radial direction and decreases near the 
black hole (see Fig. 6 in Paper 1). Thus, we conclude 
that the simulated flows do not perfectly agree with the 
slim-disk model with regard to the whole structure of the 
flow. 

iHeinzeller fc Duschll (|2007f ) have investigated local 
maximum values of the accretion rate in the supercriti- 
cal disk accretion flows. In their paper they focused only 
on the force balance in the vertical and radial directions 
around the equatorial plane. Multi-dimensional effects 
were not taken into consideration. They revealed that 
the vertical radiative force limits the maximum accre- 
tion rate at the inner disk region, leading to a decrease 
of the accretion rate with a decrease of the radius. Their 
results imply that the disk loses mass via the radiatively 
driven outflows and the mass accretion rate decreases at 
the inner region. Such tendencies agree with our results. 
As shown in Figure IH our simulations show that radia- 
tively driven outflows form above and below the disk. 
The mass accretion rate decreases with a decrease of the 
radius (see Fig. 6 in Paper I). 

3.2. Dependency of the Mass Accretion Rate 

In the present study, focusing on numerical simula- 
tions in which the mass input rate at the outer bound- 
ary is set to be Minput = 10^Fe/c^, we show that the 
radiative force is attenuated in the disk region via large 
optical thickness, which makes supercritical disk accre- 
tion possible. Such dilution of the radiative force would 
effectively operate even if the mass input rate (mass ac- 
cretion rate) varied. In fact, simulations with mass input 
rates of 3 x 10^Le/c^ and 3 x IO'^Le/c^ show that Eq/p 
is almost independent of the mass input rate, although 
the radiation energy density goes up as the mass input 
rate increases. That is, the dynamics is not sensitive to 
the precise value of Minput, as long as it largely exceeds 
the critical value. 

So far, we have studied steady accretion flows. Al- 
though a highly supercritical disk (Minput > 3 x 
10^Fe/c^) is quasi-steady, it has been revealed that 
a moderately supercritical disk [Minput = (10 — 
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Fig. 7. — Schematic picture of the supercritical disk accretion flow around a black hole. The gas motion is shown on the right: the 
high- velocity outflow {solid arrows) and the slow accretion flow {dashed arrows). The left side indicates the radiative fluxes in the comoving 
frame {black arrows) and the rest frame (white arrows). 

10^)Le/c^] is unstab le and ex hibits limit-cycle behavior 
jghsuga"'2006', '200^^;_see also IShibazaki fc Hoshil [iQTSi : 
[Sha kura & Sunvacv| |l976D . The luminosity goes up and 
d own around th e Eddington luminosity. 

lOhsugal (|2007i) has reported that the time-averaged 
mass, momentum, and kinetic energy output rates via 
the outflow, the mass accretion rate, and the disk lu- 
minosity increase as the mass input rate increases, oc 

rO.6 



<put - ^^input for a = 0.5 and (x M^^^^, 
a = 0.1. 



My-^ for 



3.3. Future Work 

As we have already mentioned in §3.1, the sum of the 
accreting matter and the matter ejected as high-velocity 



outflows is 20% of the injected mass, and 80% of the 
injected matter is ejected from the computational do- 
main as low- velocity outflows, whose velocities do not 
exceed the escape velocity at the outer boundary. Since 
such outflowing matter tends to be accelerated by the 
radiative force even at the outside of the computational 
domain, it would be blown away from the system. How- 
ever, a part of the outflowing matter might return to 
the vicinity of the black hole through the disk region, 
since the radiative force does not exceed the gravity near 
the equatorial plane. Numerical simulations with larger 
computational domains would make this point clear. 

Whereas the resulting mass accretion rate onto the 
black hole is around 10^Le/c^ in the present simulations. 
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iHeinzeller fc Duschll (|2007f ) have indicated that the mass 
accretion rate can increase up to 10^Le/c^- They have 
reported that the vertical force balance breaks down via a 
strong radiative force if the mass accretion rate exceeds 
this limit. However, even in such a case, the matter 
might accrete onto the black hole, although the strong 
radiative force would produce powerful outflows. To in- 
vestigate the maximum value of the accretion rate is an 
outstanding issue. We should perform numerical sim- 
ulations with larger computational domains, since the 
trapping region is expected to expand with the increase 
of the mass accretion rate. 

We reveal in the present paper that photons gener- 
ated deep inside the disk are effectively trapped in the 
flow, leading to supercritical disk accretion. Although 
the magnetic fields are not solved in our simulations, 
magnetic buoyancy might play an important role in the 
transportation of matter, as w e ll as photons, toward 
the disk surface (iParkeJ 1X9751 : iStella &: Rosned 1X9841 : 
iSakimoto fc CoronitillX989l ). Magnetic buoyancy might 
lead to photon generation near the disk surface if the 
magnetic fields rise quickly without dissipation deep in- 
side the disk. Thus, magnetic buoyancy would dilute 
the photon-trapping effect. A photon bubble instability, 
which is induced in the magnetized, radiation-pressure- 
domi nated region, might also suppress t he photon trap- 
ping (|Begelmanll2002l : lTurner et al.ll2005l ) . In these cases 
the enhanced radiative force would more effectively ac- 
celerate the matter around the disk surface, working to 
decrease the mass accretion rate. However, the magnetic 



Abramowicz, M. A., Czerny, B., Lasota, J. P., & Szuszkicwicz, E. 

1988, ApJ, 332, 646 
Balbus, S. A. 2003, ARA&A, 41, 555 
Begelman, M. C. 1978, MNRAS, 184, 53 
Begelman, M. C. 2002, ApJ, 568, L97 
Boiler, T. 2004, PThPS, 155, 217 
Burger, H. L. & Katz, J. I. 1980, ApJ, 236, 921 
Burger, H. L. & Katz, J. I. 1983, ApJ, 265, 393 
Ebisawa, K., Zycki, P., Kubota, A., Mizuno, T., & Watarai, K. 

2003, ApJ, 597, 780 
Eggum, G. E., Coroniti, F. V., & Katz, J. I. 1987, ApJ, 323, 634 
Fabbiano, G. 1989, ARA&A, 27, 87 

Hawley, J. P., Balbus, S. A. & Stone, J. M. 2001, ApJ, 554, L49 
Heinzeller, D., & Duschl, W. J. 2007, MNRAS, 374, 1146 
Hirose, S., Krolik, J. H., & Stone, J. M. 2006, ApJ, 640, 901 
Houck, J. C. & Chevalier, R. A. 1991, ApJ, 376, 234 
Icke, V. 1980, AJ, 85, 329 

Jones, B. C. & Raine, D. J. 1979, A&A, 76, 179 

Kato, S., Fukue, J., & Mineshige, S. 1998, Black-Hole Accretion 

Disks (Kyoto: Kyoto Univ. Press) 
Kawaguchi, T. 2003, ApJ, 593, 69 
Kley, W. 1989, A&A, 222, 141 
Kley, W. & Lin, D. N. C. 1999, ApJ, 518, 833 
Levermore, C. D. & Pomraning, G. C. 1981, ApJ, 248, L321 
Makishima, K., et al. 2000, ApJ, 535, 632 

Machida, M., Matsumoto, R., & Mineshige, S. 2001, PASJ, 53, LI 
Maraschi, L., Reina, C, & Treves, A. 1976, ApJ, 206, 295 
Meier, D. L. 1979, ApJ, 233, 664 



fields might prevent such acceleration if they strongly 
tie the matter near the disk surface with the disk mat- 
ter. In the disk region the matter might easily accrete 
toward the black hole, since the radiation energy den- 
sity decreases. Global radiation-magnetohydrodynamic 
(RMHD) simulations would make these points clear. Lo- 
cal RMHD simulations of accr etion flows have been per - 
formcd bv iTurner erall (|2003[ ) and lHirose et al.l (|2006D . 
In addition , it is thought tha t disk viscosity has ma gnetic 
origins yjawlev et al. '2 0011: iMachida et al.l 120011: for a 
review see Balbus 2003). Hence, we should stress again 
that RMHD simulations are very important to more real- 
istically investigate viscous accretion flows, although an 
a-viscosity model is employed in the present study. 
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